Negative fluctuation-dissipation ratios in the backgammon model 
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We analyze fluctuation-dissipation relations in the Backgammon model: a system 
that displays glassy behavior at zero temperature due to the existence of entropy 
barriers. We study local and global fluctuation relations for the different observ- 
ables in the model. For the case of a global perturbation we find a unique negative 
fluctuation-dissipation ratio that is independent of the observable and which diverges 
linearly with the waiting time. This result suggests that a negative effective temper- 
ature can be observed in glassy systems even in the absence of thermally activated 
processes. 

PACS numbers: 



I. INTRODUCTION 

Understanding nonequilibrium systems remains one of the major open problems in mod- 
ern physics. In the last years many theoretical and experimental studies have focused on 
the extension of the concept of temperature to the nonequilibrium regime jl| . 

Glassy systems are adequate for testing nonequilibrium generalizations of thermodynamic 
concepts. Glassy materials display extremely slow dynamics as they approach the amorphous 
solid phase from the liquid phase [2]. Below the glass transition temperature, relaxation 
times become huge and time-translational invariance (TTl) is lost meaning that two-time 



correlation and response functions strongly depend on the time elapsed since the system 
was prepared in the nonequilibrium state. At equilibrium, linear response and correlation 
functions are related by the fluctuation-dissipation theorem (FDT) j^. Although FDT does 
not hold under nonequilibrium conditions, it can be generalized by defining an effective 
temperature J]: 

Tes{t,tw) = t:,i.\ \ ^ > 5 (1) 

where C{t,tyj) is a generic two-time correlation function and R{t,tm) is the corresponding 
response of the system to an external perturbation applied at a given previous time t^. At 
equilibrium T^s is just the bath temperature. But, what is the true physical meaning of the 
nonequilibrium Tefj(t, t^)? Can it be used to characterize the nonequilibrium relaxation? 
Is it a well defined parameter from a thermometric point of view? In the last years many 
studies have tried to answer these questions from both empirical and theoretical perspectives. 
However, there are still several debated issues (for a review see ref.jsj] and references therein). 
The effective temperature is often expressed in terms of the so-called fluctuation-dissipation 
ratio (FDR): 

X{t,t^)= f t>t^ . (2) 

X{t, tw) = 1 for systems at equilibrium. In general, the asymptotic value of the FDR does 
depend, not only on the nature of the system but also on the type of perturbation applied . 
A property that a physically meaningful effective temperature Tefj(t, t^) should satisfy is 
its independence on the type of observable used to define the correlation and conjugated 
response functions in the limit t ^ t^. A standard way to account for possible differences 
is to calculate or measure X{t^tw) for different observables to evaluate such independence. 

In order to analyze the applicability and generality of the concept of effective temper- 
ature, a variety of exactly solvable models with glassy dynamics have been studied in the 
last years. A remarkable aspect of glassy systems is the appearance of negative effective 
temperatures under some conditions. This seems to contradict our intuition and to pre- 
clude a possible thermometric interpretation of the effective temperature. Recent studies 
on kinetically constrained models reveal negative FDRs % |8|] which have been interpreted 
as due to activated effects in the dynamics of such class of models. Negative FDRs seem to 



be unrelated to any thermodynamic interpretation of effective temperatures. However, from 
a theoretical point of view, nothing prevents that they could be generally found in glassy 
materials. 

In the present paper we study FDRs in the context of the Backgammon model (BG) j^. 
The BG at low temperatures presents the typical behavior of the nonequilibrium relaxation 
of structural glasses: extremely slow relaxation, time-dependent hysteresis effects, activated 
increase of the relaxation time and aging. The most interesting feature of the BG is the fact 
that glassy behavior is only due to the emergence of entropic barriers rather than energy 
barriers. 

We have found observable-independent negative FDRs in the BG due to the entropic 
barriers present at low temperatures. We conclude that the negativeness of these FDRs 
is a consequence of the dynamic coupling between the external field and the energy of the 
system. Interestingly, we also have found how these negative FDRs scale with the waiting 
time. 

The paper is organized as follows: in section [III we briefly review the BG; in section IIIII 
we present the exact analytical expressions for the correlations and responses of a set of 
correlations and conjugated responses in the model; in section [IV] we present both numerical 
and analytical results; finally, in section |V] we discuss the results. Technical aspects are left 
to the Appendixes. 



II. THE MODEL 



The Backgammon Model (BG) is a mean-field model for a glass without energy barriers. 



The model was introduced in M and has been extensively studied in lid . 

"in n 

151 . Ilq . Similarly to the case of kinetically constrained models p/7[, the statics of this 
model is very simple and does not show any phase transition at finite temperatures. The 



n 



11 



12 



BG belongs to the more general 



original Ehrenfest model 



18 



class of models called urn models which are based on the 
19| and consist of a set of M boxes (" urns" ) among which we 
can distribute particles. In these models there is no local kinetic constraint but there 
exists a conservation law, the total number of particles, that acts as a global constraint 
which induces a condensation transition. For a review of urn models and their extensions 
see Refs. 



20 



211] and references therein. 



Consider distinguishable particles which can occupy M different boxes. Let us denote 
the density (number of particles per box) by p = ^ . The Hamiltonian in the Backgammon 
model is defined by: 

M 



^ = - E ^"-0 , (3) 



r=l 



where rir is the occupation number of the box r = 1, M . The conservation of the number 
of particles gives a global constraint: 



M 



Y,nr = N . (4) 



r=l 

Eq-® shows that energy is simply given by the number of empty boxes (with negative 
sign). The system at very low temperatures tends to empty as many boxes as possible by 
accumulating all particles in a small fraction of boxes. We define the occupation probabilities 
as follows, 



which is the probability of finding one box occupied by k particles. In the canonical ensemble 
the statics can be easily solved [ol, [l^, giving the following relation for the occupation 
probabilities: 

^ ^ ^'"'eXp(/3M ^ /gN 

k\ exp(z) 

where z is the fugacity and (3 is the inverse of the temperature T. These quantities are 
related by the condition, 

p{eP-l) = {z-p)e , (7) 

expressing the fact that the density is fixed to p. This condition, in the microcanonical for- 
mulation, is equivalent to the saddle point condition in the integral solution of the partition 
function. In the grandcanonical formulation this closure condition is easier to obtain by 
means of the equation of state. The occupation probabilities are the main observables in 
the system and verify the relation YlT=o^k = 1. In particular the energy is simply given by 
U = -Po. 

In the original formulation the model was studied under Metropolis dynamics where at 
each time step a particle is chosen at random and a destination box selected. The move is 



accepted with probability one if the energy either decreases or does not change, and with 
probabihty exp(— /5) otherwise (see figure [T]). Note that the energy can only increase by 
one unit at each time step. The original geometry is mean-field, so the destination box 
is chosen at random with uniform probability among all boxes. In this case, a complete 
analytical study can be done and a hierarchy of dynamical equations can be obtained for 
the occupation probabilities |10 |. 
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FIG. 1: (Color online) Schematic representation of the dynamics of the model. At each time step 
a particle is chosen at random and a destination box selected with a uniform probability among 
all boxes. In the original formulation, the system was studied under Metropolis dynamics. 



It has been shown that the dynamics is highly non trivial at very low temperatures where 
a dramatic slowing down of the relaxational kinetics takes place. The origin of this slowing 
down can be qualitatively understood. Suppose that the system starts from a configuration 
of high energy and the temperature is set to zero. The system will then evolve without 
accepting changes which increase the energy of the system. As time goes on, the system 
evolves towards the ground state of the system where all boxes are empty and all particles 
have condensed into a single box (figure [2]). During the relaxation process more and more 
boxes are progressively emptied; this means that the few boxes which contain particles have 
more and more particles because the number of particles is a conserved quantity. Then, the 
time needed to empty an additional box increases with time. The final result is that the 
energy very slowly converges to the ground state value. At very low temperatures it can be 
shown that the characteristic equilibration time is given by: 

, exp/3 

r = teq ^ (8) 



which diverges at zero temperature. The Arrhenius dependence is remarkable if we note 
that only entropy barriers (but not energy barriers) are present in the model. 



T = 




FIG. 2: (Color online) At zero temperature only movements between occupied boxes are accepted. 
As time goes on, only a small fraction of boxes contain particles and the time needed to empty an 
additional box increases rapidly as the number of occupied boxes decreases. 



The BG has been used as a playground model where new concepts of nonequilibrium 
thermodynamics can be tested. The fact that the dynamics is glassy and can be exactly 
solved has inspired several works that have investigated extensions of FDT to the nonequi- 
librium regime (e.g. the disordered model studied in 22|). In the present work we solve the 
BG for any general Markovian dynamics and study the existence of negative FDRs. 



III. CORRELATIONS AND RESPONSES IN THE BG 

Let us generalize the BG by adding an external field to the Hamiltonian of the model ([3]). 
The external field is introduced in order to compute the effective temperature, Eq. ([T]), in 
the nonequilibrium regime. This external field can be a local quantity (i.e. an external field 
acting on a single box) or a global one (i.e. an extensive field acting on all boxes) leading 
to different definitions of the FDRs. 

Previous studies of the nonequilibrium dynamics of the BG, such as the studies carried 



out in 



15|, have suggested that the effective temperature depends on the observable. In the 



studies of Ref. 



151 ]. the Hamiltonian was perturbed by a /oca/ external field. Recently, it has 



been shown that local FDRs can lead to inconsistent results if finite-N corrections are not 
properly taken into account 23|], pointing out the convenience of computing global FDRs. 



In order to give a complete picture of the system, throughout this paper we will compute 
both local and global FDRs by considering local and global external fields. 



A. Local external field 



Let us consider an external field h acting on a single box (e.g. box number one) that is 
coupled to this box only when it contains one particle: 

N 

H = Y^ (-5n.,0) - h5^n,,l . (9) 
r=l 

Note that this subextensive perturbation does not affect the values of the occupation 
probabilities Pkif) = < J2r^nr,k > which in equilibrium are still given by eq.®- As can 
be deduced from (9) we set, without loss of generality, the density of the system p = 1. 
However, note that all the results obtained throughout the paper are valid independently on 
the value of the density p = N/M whenever p is finite in the N ^ oo limit. In Appendix 1X1 a 
complete derivation of the dynamical equations for the probability densities of the perturbed 
box is carried out (see eq. OAip ). In what follows we will focus on the dynamical evolution 
of two-time quantities: local correlations and local response functions. Local correlation 
functions are defined as: 

Ci"% t^) = ^<Yl ^n.it),kSnAt^),l > , (10) 
r 

where the sum in (ITOl) runs over all boxes and counts the fraction of boxes that contain k 
particles at time t provided that these boxes contained one particle at previous time ty^. The 
brackets denote an average over dynamical trajectories of the system and over the initial 
conditions. The dynamical equations for these local correlations are derived in Appendix |X] 
leading to (see eq. lASp : 



' ' = w{o)[-kci"^ + {k + - + ] + 

+ (1^(0) - W{-mP^{&r - Ci-i) + (4,1 - 4,o)(C^(l - Po) + C^Pi)] + 

+ iW{0) - WimPoikCi"' -{k + + (4,0 - 4,i)(C(l - Pi) + C^Po)] . 

(11) 

This expression is valid for any Markovian dynamics. W{AE) denotes the transition 
probabihty between two states with energy difference AE. From now on, we consider that 
the dynamics obeys local detailed balance in order to ensure the convergence toward equi- 
librium. 

Similarly, we can compute the dynamical equations for the local response function defined 
as the variation of the occupation probabilities for the perturbed box when the impulse field 
is applied at t^: 

Ri"'it,t^) = (^SM] . (12) 



Again, the details about the derivation can be found in the Appendix |Al The result 
(eq.dMD) is: 



dRi^'^{t, tu 
di 



= W^(0)[-fcPf +ik + l)R'°l, - Pf + Pi-,] 
+ (W(0) - iy(-l))[Pi(Pf - Pif^i) + - 4,o)(/2^(l - Po) + Ro'Pi)] 
+ {W{Q) - iy(l))[Po(A;Pf - (fc + l)R':i,) + (5fc,o - 4,i)(P^(l - Pi) + P'rPo)] 
+ 5{t-t^)S'°'\<P,>] , (13) 

where the (5-term S^°'^[< >] is given in eq. flASp . Equations ffTTl) and (fTHjl are the necessary 
ingredients for computing nonequilibrium effective temperatures. 

From ffTTl) and ffT^ . we can check that FDT is verified at equilibrium. Indeed, at equilib- 
rium the correlations and responses become functions of the difference of times, i.e Cj!"^{t—tu,) 
and R["%t-t^) so we recover time-translational invariance. Moreover, as we can see from the 
form of the dynamical equations for the autocorrelations (eq. lfTTl) ) and responses (eq.( fT3l) ). 
at equilibrium the FDT is verified at all times provided that the initial condition for the 



responses (the function S''°'^[< >]) corresponds to the value of the derivative of the ap- 
propriate correlation at equal times. In this case, the correlation functions for a general 
observable are given by: 

= Pi(t^)4,i . (14) 
Therefore, the initial value {t = t^) for the derivative of the correlation functions are: 



' ^ ' ' - Pi{t^) + 5fc,o + 4,2)] + 



dt 



+ Pi(t^) [Wio) - wi-1)] (Pi(4,i + 4,2)) + 



+ Pi(t^) [WiO) - Wi-l)] (1 - Po)(4,i + M . (15) 
Using eg. (lASp it is easy to check that, in equilibrium, FDT is verified. 

T=-—^—^ -. (16) 

B. Global external field 



As we have explained before, local computations can 



ead to erroneous conclusions if finite- 



corrections are not properly taken into account [23|]. In such cases it is easier to carry 
out an analysis of FDRs for global observables. Here we shall consider the corresponding 
extensive perturbation of an external field coupled to the set of boxes which contain just 
one particle (i.e coupled to the observable Pi). The Hamiltonian reads: 

N 

H=-J2{Sn^,0 + h6n^,,) . (17) 
r=l 

Now, as the perturbation is extensive, all the equilibrium occupation probabilities are 
modified in the presence of the external field h: 

p _ z''-^ exp{f3Sk,o - f3hSk,i) , . 

' k\{e- + e-f'^ - 1) ■ ^ 

We proceed following the same steps as in the local case. In Appendix [B] we have com- 
puted the dynamical equations for the occupation probabilities, eq. llBlll . and from these 
equations we derive the dynamical evolution for the global correlation and response func- 
tions. 



Due to the fact that the perturbation is extensive we consider the connected correlation 
functions: 



Clit,t^) =< -ik{t)-il{tn,) > , (19) 

where 

7kit) = ^J2^-r,k-Pkit) (20) 

r 

are the deviations of the instantaneous occupation variables from their average value at a 
given time. The dynamical evolution for the global correlation functions is given by eq. (lB6p : 



%; = wm-kci + {k + i)cu, - ci + + 

+ (W^(O) - W^(-l))[Cf (4,1 - 4,0 + Pk- Pk^i) + Pi{Cl - + 

+ (W^(O) - W{l))[Cf,{kPk -{k + l)Pfe+i + 4,0 - 4,i) + PoikCl -{k + l)C,\i)]. 

(21) 

Again, these equations are valid for any Markovian dynamics. The global response function 
is the response of the occupation probabilities to the extensive perturbation coupled to the 
observable Pi: 

sPk{ty\ 



The result for the dynamical evolution is given in Appendix [Bl eq. (lB8l) . and it reads: 



= wm-kRi + {k+ mu - K + RU] + 

+ (W^(O) - W^(-i))[P?(4,i - 4,0 + Pk- Pk^i) + Pi{Ri - Ri.i)] + 

+ (iy(0) - W{l))[Rl{kPk -{k + l)Pfc+i + 4,0 - 4,i) + Po{kRi -{k + l)Pf+i)] + 
+ 5{t-t^)S^[<Pu>] , (23) 

where we have introduced the function 5'^[< P^ >] which depends only on one time and 
gives the initial value for the responses. The exact form of S^[< Pk >] is given in eq. ( 1B9I) . 

Again, we check that in equilibrium FDT is verified. Indeed, at equal times the global 
correlations are given by: 



Clit^, t^) = -PMiPiitu,) - 4,1) • (24) 

Obviously, in equilibrium the correlations at equal times do not depend on time. Inserting 
this initial value in the equations for the correlations and by considering the equilibrium case, 
it is easy to prove FDT for all k values of the observables C^, Rl: 

dCl{t-t^) 

IV. RESULTS 



In this section we analyze the nonequilibrium behavior of the correlations and responses 
at zero temperature both for local and global observables. 

The interesting glassy behavior in the BG model occurs in the zero temperature limit 
where entropy barriers govern the relaxational dynamics of the model. In what follows we 
shall consider heat-bath dynamics at zero temperature, both for the local and the global 
variables. This choice is motivated by the known fact that in the Metropolis algorithm 
there is a discontinuity of the derivative of the transition rates for AE = 0. As a result, 



the definition of the response functions become ambiguous, see Ref. 
drawback by employing heat bath dynamics. 



15|. We circumvent this 



A. Local two-time quantities 

From the numerical integration of equations ( fTTl) and ( fT3l) we can analyze the nonequi- 
librium behavior of the local correlations and response functions. From now on, all the 
numerical results shown are obtained using heat bath dynamics at zero temperature. 



1. Correlations and responses 

In figure [3] we plot the normalized local correlation &°^{t,tyj) = '"p^^^^™^ at zero tempera- 
ture. We can clearly see the aging effects in the local correlation function: as tyj increases the 
autocorrelation function develops a plateau showing two characteristic and well separated 
time-scales: the first timescale corresponds to the initial relaxation of the system (usually 



called /3-relaxation) which does not depend much on t^; the second one is larger, increases 
with and corresponds to the late decay of the correlation function, usually known as 
a-relaxation. The existence of these two time-scales is a typical signature of the glassy 
relaxation of structural glasses. 

In the inset of figure [3] we plot the local normalized correlation function C'{°'^(t, t^) mul- 
tiplied by in order to collapse all curves on the same plateau. It is clear that the system 
displays simple aging, i.e the scaling t/tyj is well satisfied. 
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FIG. 3: (Color online) The evolution for the normalized local correlation function C[°'^{t, t^) (di- 
mensionless) for different t^'s. In the inset the scaling of C'{'"^(t,t^) corresponding to simple aging 
is displayed. The time is measured in Monte Carlo sweeps. 



Regarding response functions, they show some peculiarities: on the one hand, the initial 
value for the response functions (given by the function S''°'^[< >]) is proportional to /3, 
giving a divergence at zero temperature (a known common feature of kinetically constrained 
models Q]); on the other hand, the response function B!'°^{t,tw) is non-monotonic (for 
fixed when t is varied) and becomes negative for long enough times. In figure H] we plot both 
the local, B}°'^, and the global, response functions for the observable Pi (see below). 
Both responses show a non-monotonic behavior and become negative for long times. 

The non-monotonicity of the response function can be easily understood. The external 



field is coupled to Pi, therefore the system tends to increase the population of boxes with 
one particle; however, because boxes with one particle are bottlenecks for the relaxation 
of the energy, a transient increase in their number at t^, induces a faster relaxation of the 
energy at later times. Because the natural evolution of the system tends to decrease Pi 
when decreasing the energy, a transient increase of Pi at ty^ causes a net decrease of the 
same quantity at later times when energy relaxation becomes faster. 

In order to facilitate the readings of the effective temperature in our plots we introduce 
the function G[°%t, t^) as: 



Gf(t,t^)=T 



\{Rr{t,t^))\ 



(26) 



which is the normalized absolute value of the response B}^^{t,t^) muhiplied by T. Due to 
the change of sign of B!'"'^, we have taken the absolute value in order to plot the relaxation 
in a log-log scale. 

In figure E] we plot G^°'^{t,tw) for different values of t^. As can be inferred from figure El 
the dip at short times corresponds to the change in sign of the response. Looking at this 
logarithmic plot, the response shows again the two characteristic relaxation time-scales of 
glasses. In the inset of figure [5] we can see that the response function also displays simple 
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FIG. 4: (Color online) Time dependence of the global and local dimensionless response functions 
multiplied by T {TRf and TR'f'^) at tyj = 1000. The time is measured in Monte Carlo sweeps. 
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FIG. 5: (Color online) The evolution of the dimensionless quantity, G''f'^{t,tw), (|26|) at different 
values of ty^. In the inset we show the simple aging scaling for this quantity. The dips observed in 
qIoc infiicate changes of sign in Rlf^. The time is measured in Monte Carlo sweeps. 

aging with scaling t/tw as the leading term. 

It is worth noting that this simple aging relaxation in the a-regime can also be seen in 
the correlations and responses for the other observable quantities of the model, i.e in the 
dynamical behavior of C^^'^(t,ty^) and E}^'^{t^tw) for a generic k (data not shown). 

2. Nonequilibrium effective temperatures 
We now define a set of effective temperatures from the nonequilibrium definition, eq.([l]): 

In order for (T^'j^^^)^ (t, t^) to share some of the properties of a thermometric temperature 
it should not asymptotically depend on the integer k (for a fixed and in the large t limit). 
In figure E] we plot the ratio between the absolute value of the effective temperature and the 




FIG. 6: (Color online) The absolute value of the effective temperature divided by T for the observ- 
able Pi (/c = 1) as a function of time for different t^'s. The (up and down) oriented spikes indicate 
changes in the sign of (T'^g'^)^. Note that, for a given t^, the effective temperature changes sign 
twice. Therefore, we can distinguish three regions depending of the sign of the local FDR. The 
time is measured in Monte Carlo sweeps. 



physical one (in the limit T — 0), which corresponds to the inverse of the local FDR ^ 
defined as: 



We can clearly see that the effective temperature shows two different behaviors depending 
on the time-scales considered. For t t^, the value of the FDR converges to 1 as t^, increases. 
This is a typical feature of glasses: the first /^-relaxation is an equilibrium process which 
implies that the effective temperature is just the physical one. This is true in the asymptotic 
limit oo. It can be shown that it converges to 1 in a logarithmic way as was found in 

the analysis of Ref . [isl . 

From the integration of the dynamical equations, we obtain the asymptotic value of the 
ratio (Tgff"^) ^ (t, tu,)/T, which is positive because for long enough times both the local response 
and the derivative of the local correlation become negative. This asymptotic value tends to 
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FIG. 7: (Color online) The absolute value of the effective temperature divided by T = 10000) 
as a function of time for different observables k. In the inset we zoom the boxed part of the figure 
in order to emphasize the observable dependence of the effective temperature. The (up and down) 
oriented spikes indicate changes in the sign of {Tl,g'^)j^. The time is measured in Monte Carlo 
sweeps. 

zero in the limit ^ oo. In addition, for a given waiting time the effective temperature is 
proportional to the bath temperature. 

Looking at figure [3, where we plot the effective temperature at = 10000 for different 
observables, we can clearly see that the effective temperature depends on the observable 
under scrutinity. Consequently, it seems clear that from a local point of view we cannot 
define a unique effective temperature by using the FDR. 

B. Global two-time quantities 

In the preceding section we have shown that a unique effective temperature cannot be 
defined by the FDR from a local perturbation. As we have mentioned before, this is an 
expected result consistent with previous 

time dependence of the global correlation and the global response functions and we will show 
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FIG. 8: (Color online) The global (connected) correlation function normalized for different t^'s. 
In the inset we plot the simple aging scaling for Cf{t,tw)- These correlation functions are dimen- 
sionless and the time is measured in Monte Carlo sweeps. 

that a unique negative effective temperature can be defined from the global FDRs. 

1. Correlations and responses 

We study the connected correlation functions for heat bath dynamics of the BG at zero 
temperature. In FigIH]we plot the normalized correlation function, Ci{t,t^) = -^jj-y-, for 
different values of t^. Similarly to the local case, we can clearly distinguish two characteristic 
time-scales in the system, the /3-relaxation and the a-relaxation. Note that, as increases, 
the plateau value of the correlation decreases and in the limit t„, —>■ oo the plateau value 
converges to zero. In the inset of FiglH] we have multiplied this normalized correlation by 
tuj. As for the local case, the global correlation displays simple aging. 

Again, in order to analyze the relaxation of the global response function we have defined 
the normalized response function Gl(t,t^) 
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FIG. 9: (Color online) The global scaled response function for Pi for different f«,'s. The dips 
observed in (dimensionless) indicate changes in the sign of Rf. The time is measured in Monte 
Carlo sweeps. 



G?(t,U = T^^^^^^, (29) 

motivated by the fact that the response function Rf{t, ty^) becomes negative for long times as 
it is shown in figure HI This is again consequence of the fact that the natural evolution of the 
system tends to diminish Pi(t^) in opposition to the action of the external field. Moreover, 
the global response function is proportional to the bath temperature, which diverges at zero 
temperature. In figure [9] we plot the two time-scales relaxation of G\{t^t^). In the inset of 
figure [9] we show the simple aging scaling of the function Gf{t,tw)- Again, the dip of the 
curves at short times corresponds to the time when the response changes its sign. 

The global correlations and responses for the rest of the observables in the model, Cl{t, t^) 
and Gl{t,tw), also display simple aging (data not shown). It is worth mentioning that the 
ratio between and for /c > 1 is of the same order of magnitude as the one corresponding 




FIG. 10: (Color online) The absolute value of the global effective temperature function normalized 
by T for Pi for different values of t^. The up oriented spikes indicate changes in the sign of (Tl^gj^ 
which is negative for long times. The time is measured in Monte Carlo sweeps. 



2. Nonequilibrium effective temperatures 



As we have done for local observables, from the FDR we can define the effective temper- 
atures: 

ra.(^.u = ^^- (30) 

In figure [To] we plot the absolute value of the effective temperature (Tjj)-^ (t,^^) divided 
by T for different values of tw This quantity is related to the global FDR Xfit^) for Pi as: 

As in the local case, the value of °^ \^ in the limit t ^ tends to 1 showing that the 
first /3-regime corresponds to an equilibrium relaxation process (i.e Xf = 1). We can also 
see that, in contrast with the local case, this global effective temperature remains constant 
throughout the a-regime for any finite t^. 
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FIG. 11: (Color online) The global effective temperature function normalized by T for tw = 10000 
and for different observables. The (up and down) oriented spikes indicate changes in the sign of 



The time is measured in Monte Carlo sweeps. 

A very important aspect of the global effective temperature is the fact that, for a given 
value of the waiting time this effective temperature does not depend on the observable as 
can be seen in figure [TT] where we have plotted the absolute value of the inverse of the global 
FDRs (eq. flHT]) ) at t^, = 10000 for different observables. It is clear that the asymptotic value 
of the FDRs at finite does not depend on the observable. 

Moreover, from the results of figure [TU] we can see that for large waiting times t^, the 
inverse of the FDR scales as the inverse of t^: 



Again, the minus sign in eq. flH^ is a consequence of the non-monotonicity of the response 
functions. Finally, it is worth mentioning that we have checked that all the results obtained 
at zero temperature throughout this paper remain valid at finite but very low temperatures. 
The analysis at finite small temperatures does not give new insights into the nonequilibrium 
behavior of the system as all dynamical quantities smoothly converge to their T = limit. 



some values of k the effective temperature changes sign more than twice. 




(32) 



T 



C. Asymptotic analysis 



In the preceding section we have obtained a negative FDR independent of the observable 
that displays simple scaling of the type t/tw This result can be easily understood by ana- 
lyzing the asymptotic nonequilibrium relaxation of the model. The equilibrium probabilities 
in the presence of an external field are given by: 

2;'=-^exp(/5(4,o - ^4,1) /oQ\ 

• ^^^^ 

In the long time asymptotic regime, the multiplier z{t) is a function of time which grows 



as 



11 



12|: 



z{t) ^ Int + ln(lnt) . (34) 

With the global perturbation considered along the paper we can compute the global 
susceptibility Xi by assuming local equilibrium using (fT8l) with k = 1: 

T,; = TUn>^M^^ = l . (35) 
From (IHll) the global susceptibility associated to the observable Pi{t) decays as: 



= 7h^ ' ^^^^ 
where r = t—ty^. The asymptotic decay of R\{t) is given by the derivative of Xi^t) multiplied 
by the temperature: 



Now, by using the equation (l23l) at zero temperature: 



dr 

we obtain the asymptotic decay of R^: 



^^0 ^ ^9 _ p^j^, _ p^^s ^ (38) 



In the right column of figure [T2] we numerically confirm the scalings (1371) and fl39p for 
different values of t^. Due to the fact that the dynamical equations for the correlations are 
formally identical to those for the response functions one finds: 

C^f(-.«»)-^^ + f(-J^) , (40) 

where the dependence on t^^ has been inferred from the decay of the global correlations at 
equal times (eq. flM|) ). 

These scalings are again confirmed numerically and are shown in the left column of figure 
[T21 With these scalings we recover the following FDRs: 

Xo^(tJ = (t J ^ -t^ , (41) 
in agreement with our numerical findings. A similar analysis can be done for k > 1. 

V. CONCLUSIONS 

In this paper we have solved the relaxation of the correlations and response functions in 
the BG model for a general dynamic rule (provided that it satisfies local detailed balance). 
We have studied (by means of numerical integration and analytic asymptotic expansions) 
the behavior of effective temperatures and FDRs in the glassy regime. 

We have found that both the correlation and the response functions show two character- 
istic time-scales: a first /5-relaxation for short times characterized by an equilibrium FDR, 
X(t,tw) = 1, and a second a-relaxation at long times with a non-trivial value of the FDR. 
This is a very common feature of structural glasses and other glassy systems. Moreover, we 
have found that both the correlations and responses display simple aging. 

In this paper we have analyzed the resulting FDRs obtained from local and global per- 
turbations. The interesting conclusion is that the local FDRs depend on both t and 
while the global FDRs only depends on t^. Moreover, global FDRs are independent of the 
observable in contrast with the local ones. 
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FIG. 12: (Color online) Left column: scaling of the global correlations Cq and Cf divided by In(t^). 
Right column: scaling of the global correlations Rq and Rf multiplied by T. The continuous lines 
are the assymptotic scalings ([10|) and ([39]) . ([37]) while the discontinuous lines correspond, from top 
to bottom, to tw= 100, 1000, 10000 respectively. All quantities are dimensionless and the time is 
measured in Monte Carlo sweeps. 



More interesting is the fact that this observable-independent value of the global FDR is 
negative and diverges with the waiting time as: 

• (42) 

This result points in the same direction as recent studies on kinetically constrained models 
71, |8| which also found negative FDRs. In these studies, the negative character of the FRDs 
was associated to activation effects in the dynamics. In the present case, we have found 
negative FDRs in the sole presence of entropic barriers for the BG. 

It is worth emphasizing that negative FDRs are related to non-monotonic response func- 
tions. In the glassy literature, non-monotonic responses are associated with non-neutral 
observable quantities The non- neutrality property of these observables emerges as a 



consequence of the dynamic coupling between the external field and the energy of the sys- 
tem leading to negative effective temperatures. 

In this paper we showed that an observable-independent FDR can be properly defined 
by studying global observables. However, we found a unique negative FDR due to the non- 
neutrality of the observables under scrutinity. Therefore, the neutrality of an observable 
seems to be a key aspect in order to define nonequilibrium effective temperatures. 

How much current results would change if the perturbation h6n,k acts along an arbitrary 
direction k > 17 We do not expect big qualitative changes in our results depending on the 
"orientation" of the field provided that k is finite (and k « N). Arbitrary values of k will 
result in a bottleneck effect similar to that observed in the current study for k = 1. However, 
for k/N finite the bottleneck effect will be substantially different because the energy will 
not be able to reach the asymptotic low energy regime E ^ —1 + . 

Finally, it would be extremely helpful to find a microcanonical derivation or a phenomeno- 
logical argument for reproducing the asymptotic behavior of the effective temperature when 
perturbing along arbitrary observables Pk- This could be done either by a closure of the 
dynamical equations by using a partial equilibration hypothesis, or by exact computation 
of the appropriate configurational entropy in the off-equilibrium regime. Such arguments 
would greatly facilitate the computation of effective temperatures without having to solve 
the full set of dynamical equations for correlations and responses. 
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APPENDIX A: LOCAL DYNAMICAL EQUATIONS 

In the present analysis we consider a general dynamics with just one restriction: it must 
obey local detailed balance. This restriction ensures that the system converges to its equilib- 
rium state. In fact, we will see that this is the necessary condition for FDT to be obeyed at 



equilibrium. From now on, the transition probabilities will be expressed by W{AE) where 
AE is the energy difference between the final and the initial states. 



a. One-time quantities. 



The dynamic equations for the occupation probabilities can be computed in the same 
way as have been obtained for the Monte-Carlo dynamics of this model (see Ref. HI] for 
details). The general result is: 



A = wm-kp^ + ik + - H + Pk-i] + 

+ {W{0) - W{-1 + h))[Pl{l - Po){6k,i - 6k,o) + - 4,2)] + 

+ {W{0) - W{h))[PlPo{6k,i - dk,o) + - Pim,i - 4,2)] + 

+ iW{0) - Wi-h))[2P^{l - Po){6k,2 - 6k,i) + P^Pi{6k,o - 6k,i)] + 

+ {W{0) - W{1 - h))[2P^Po{6k,2 - 4,1) + 4(1 - Pi){Sk,o - 4,1)] + 

+ {WiO) - WimPoikP^ -{k + l)Pl_^_,) + PlPo{Sk,o - Sk,i) + 2P^PoiSk,i - 6k,2)] + 

+ {W{0) - W{-l))[Pi{Pt - P^) + P^Pi{6k,i - Sk,o) + PlPi{6k,2 - . 

(Al) 



The quantities Pk{t) are the occupation probabilities, while the quantities Pl{t) are the 
average occupation probabilities restricted to the box number one, which is the box af- 
fected by the external field. As a particular case we can get the dynamic evolution for the 
occupation probabilities for any box at zero field: 



dt 



= W{Q)[-kPu + {k + l)Pfe+i - Pu + Pfe_i] + 

+ {W{Q) - W{-l))[Pi{Pk - Pk-i + 4,1 - 4,0)] + 

+ {W{0) - W{l))[Po{kPk -{k + l)Pk+i + 4,0 - 4,1] 



(A2) 



These equations cannot be solved exactly (although an analytic treatment has been done 



in the asymptotic regime 



but can be integrated numerically to give the full solution. 



More significative, these equations are the first step in order to compute the dynamical 
evolution of two-time quantities such as the autocorrelation functions and the local response 
functions. 



b. Local correlations and response functions 



As a consequence of the local character of the external field, we have to deal with the 
corresponding local response functions and the box-box autocorrelation functions. These 
autocorrelation functions are defined as follows: 

Ck''{t,t^) = — <^6n,{t),kSnritu,),l > ) (A3) 
r 

which can be expressed in terms of the following conditional probabilities Vkititw) = 
Pijirit) = k\nrit^) = 0): 

Cf(t,t^)=Pi(Uz/fc(t,U • (A4) 
Following the same strategy as in ll| the dynamic equations for these conditional prob- 



abilities give: 

^^^^iM^) = W{0)[-ku, + {k + l)z/,+i - z/, + z/,„i] + 

+ iW{0) - Wi-l))[P,{u, - z/fc„i) + {6k,i - - ^o) + ^^o^i)] + 

+ (W(0) - W{l))[Po{kuk -{k + + (4,0 - 4,1) K(l - Pi) + ^iPo)] ■ 

(A5) 

The corresponding local response functions are just the variations of the occupation 
probabilities for the perturbed box with the external field: 

From this expression and eq. (IA2p we arrive at: 



dt 

+ (W{0) - iy(-l))[Pi(i?f - Ri%) + iSk,i - 4,o)(i?'r(l - Po) + R^Pi)] 

+ {w{Q) - iy(l))[Po(fci?f - (A; + l)4+i) + (5fe,o - - Pi) + P'rPo)] 

+ 5(t-U5'°'=[<Pfc>] (A7) 

Note that formally, the dynamic evolution for the response functions is just the same 
as for the autocorrelation functions. This is a general feature and is due to the fact that 
in equilibrium FDT must be satisfied. The only difference is that in the equation for the 
responses there is a delta term which fixes the value for R^j^'^{tu],tw). This term comes 
from the first order of the Taylor expansion in the transition probabilities which depend on 
the external field h. This is not an approximation because higher order terms in Taylor's 
expansion vanish when we set the external field equal to zero. The function ^'°=[< Pk >] is 
defined as: 

S''"'[< Pk >] = Pe^W{l)[P,il - Po)(4,i - 4,0) + Pi'(4,i - 4,2)] + 
+ /3e^iy'(l)[Pi(l - Po)(4,i - 4,0) + Pi'(4,i - 4,2)] - 

- (3W{0)[2P2{1 - Po)(4,2 - 4,1) + PiPo(4,o - 4,1)] - 

- f3W'{0)[2P2{l - Po)(4,2 - 4,1)] + 

+ l3W'{l)[2P2Po{6k,2 - 4,1) + (1 - Pi)Po(4,o - 4,1)] - 

- /3W' {0)[{1 - Pi)Po(4,o - 4,1)] , (A8) 

where W denotes the derivative of the transition probability with respect to AE. Finally, we 
must stress that in this equation we have already supposed that our dynamics verifies local 
detailed balance. If W is discontinuous we should have to consider two possible response 
functions depending on the chosen value for W [l^. 



APPENDIX B: GLOBAL DYNAMICAL EQUATIONS 



As we have made in the case of a local perturbation, we consider a general dynamics 
with the only condition that it must obey detailed balance. As before, this is the unique 
ingredient we need to ensure that equilibrium is reached at long enough times. 



c. One-time quantities. 



By considering all the possible elementary moves, we get the following dynamical equa- 
tions for the occupation probabilities: 



dPk 

at 



= W{0) 
+ {W{0 
+ {W{0 
+ {W{0 
+ {W{0 
+ {W{0 
+ {W{0 



-kpk + {k + l)pk+l -Pk + Pk-l] + 

- W{h - l))[Pi(l - Pi)(4,i - 4,0) + PkPi{l - - Pk-iPi{l - 4,2)] + 

- W{2h - l))[Pi'(24,i - 4,0 - 4,2)] + 

- W{-h))[2P2{l - Po)(4,2 - 4,1) + 2PfcP2(l - 4,0) - 2P,_iP2(l - 4,1)] 

- W{1 - 2/i))[-2PoP2(24,i - 4,0 - 4,2)] + 

- W{1 - h))[Po{kPk -{k + l)Pfc+i + 4,0 - 4,1 + 2P2(24,i - 4,0 - 4,2))] 

- W{h))[P^{kPk -ik + l)Pfe+i + 4,1 - 4,2) + Pi'(-24,i + 4,0 + 4,2)] . 

fBll 



Note that now, due to the global character of the field, we have to consider only the 
occupation probabilities averaged over the whole system. As we expect, at zero field we 
recover the same equations as in the local case (which are the extension of the equations 
obtained for Monte Carlo dynamics [lol|): 



= w{0)[-kPk + {k + l)Pfc+i - Pk + Pk-i] + 

+ {w{o) - iy(-i))[Pi(4,i - 4,0 + Pk- Pk-i)] + 

+ {W{0) - W{l))[Po{kPk -{k + l)Pfc+i + 4,0 - 4,1)] . (B2) 

These equations are the first step in order to compute the dynamical equations for the 
correlation and response functions and give the evolution of all the possible observable 
physical quantities of this model. Moreover, these equations are the base of the more complex 
computations of the dynamical evolution of the two-time correlation and response functions. 



d. Global correlations and response functions 



Due to the extensive nature of the perturbation the correlation functions related with the 
responses are the connected ones. So, let us introduce the deviation of the instantaneous 
values of the occupations from their average value at each time: 

7fcW = ^5^'5n.,fc-PA.(t) , (B3) 

r 

These quantities will give us insight on the fluctuations of the occupation numbers (i.e 
the correlations). The dynamical evolution of these quantities is: 

^ = W{0)[-kjk + {k + l)7fc+i - 7fc + 7fc-i] + 

+ (W^(0) - iy(-l))[7i(4,i - 4,0 + Pk - Pk-i) + Pi{lk - Ik-i)] + 

+ (iy(0) - W{l))[^^{kPu -{k + l)Pu+i + 4,0 - 4,1) + Po{k^k -{k + l)7fc+i)]. 

(B4) 

In this equation we have considered that the quantities 7^ are of order so we have 
neglected the quadratic terms 7^7; in these equations because they vanish in the thermody- 
namic limit. The global connected correlation function will be: 

Clit.t^) =< 7fc(t)7i(iu,) > . (B5) 

The equations of motion for these correlations are easy to compute from these equations 
and give: 

^-^1^ = wm-kci + {k + i)ci^, - ci + cu\ + 

+ (W^(O) - (4,1 - 4,0 + Pfe - Pk^i) + Px{Ci - Cti)] + 

+ (H^(O) - W{X))\Cl{kPu -ik + l)Pk+i + 4,0 - 4,1) + PoikCl -ik + 

(B6) 

Now we define the global response function (which is related to the experimental sus- 
ceptibility) as the response of the probabilities to the extensive perturbation coupled to 
Pi- 



From the equations in a field and expanding to first order in h we have for the response 
functions: 



g; = wm-kRi + {k+ i)RU, - Ri + RU] + 

+ {W{Q) - W{-l))[R{{5u,i - 5m + Pk- Pk^i) + PiiRl - RLi)] + 

+ {W{0) - W{l))[R',{kPk -{k + l)Pfc+i + 4,0 - 4,1) + PoikRl -{k + l)Rl^,)] + 

+ 6{t-t^)S^[<Pk>] , (B8) 

where we have defined the function S^[< Pk >] which depends only on one time and gives 
the initial value for the responses. Similar computations as we have done for the local case 
lead to: 



S3[< Pk >] = PePW{l)[P,il - Pi)(4,i - Sk,o) + PkPiil - 4,1) + Pk^iPiil - 4,2)] + 
+ Pe^W'{l)[Pi{l - Pi)(4,i - 4,0) + PkPiil - 4,1) + Pk-iPi{l - 4,2)] + 
+ 2Pe''{W{l) + W'{l))[P^{25k,i - 4,0 - ^^,2)] - 

- PW{0)[2P2{1 - Po)fe - 4,1) + 2PfcP2(l - 4,0) - 2Pfc_iP2(l - 4,1)] - 

- /3W(0)[2P2(1 - Po)(4,2 - 4,1) + 2PfcP2(l - 4,0) - 2Pfc_iP2(l - 4,1)] + 
+ 2/3iy'(l)[-2PoP2(24,i - 4,0 - 5fc,2)] + 

+ /3W'il)[Po{kPk -{k + l)Pfc+i + 4,0 - 4,1) + 2PoP2(24,i - 4,0 - 4,2)] - 

- f3W'{0)[Pi{kPk -{k + l)Pfc+i + 4,1 - 4,2) + Pi'(-24,i + 4,0 + 4,2)]. 

(B9) 

As before W'{AE) is the derivative of the transition probability with respect to AE 
evaluated at AE = 0, 1. 
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